Short Range Interaction Effects on the Density of States of Disordered Two 
Dimensional Crystals with a half— filled band 
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The Density of electronic States (DoS) of a two-dimensional square lattice with substitutional 
impurities is calculated in the presence of short-range electron-electron interactions. In the middle 
of the energy band, the Bragg reflections off the Brillouin zone boundary are shown to lead to 
additional quantum corrections to the DoS, the sign of which is opposite to the sign of the Altshuler- 
Aronov's logarithmic correction. The resulting quantum correction to the DoS at half-filling is 
positive, i.e. the DoS increases logarithmically as the Fermi energy is approached. However, far 
from the commensurate points where the Bragg reflections are suppressed, the negative logarithmic 
corrections to the DoS survive. 
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The competition of the effects of interaction and ran- 
domness in one-dimensional (ID) and two-dimensional 
(2D) electronic systems is one of the most intriguing 
problems in low temperature physics. Impurities of arbi- 
trary concentration in a ID metal off half-filling have 
been shown to localise all electronic states of a non- 
interacting electron gas, [1]. However, the Bragg reflec- 
tions appearing at half-filling delocalise the states and 
enhance the density of states (DoS) in the middle of the 
band, an effect known as 'Dyson singularity', [2,3]. A 
ID correlated electron gas without disorder is described 
by the Luttinger liquid theory, [4] . In a half-filled band, 
charge excitations with gap (Mott insulator) may be rel- 
evant in ID systems as a consequence of commensura- 
bility, [5], whereas the localization length for coherent 
propagation of two interacting particles in a ID random 
electronic system has been shown [6] to be larger than the 
one-particle localization length, i.e., interaction leads to 
a significant derealization of the pair. 

The problem of the interplay of interaction and dis- 
order in 2D systems is as complicated as it is in ID 
systems. Indeed, all states of a 2D electron gas have 
been proved to be localized irrespective of how small the 
impurity concentration is, [7]. Weak Coulomb interac- 
tions between electrons moving diffusively in a 2D dis- 
ordered metal away from half-filling have been shown to 
increase the localization effect, [8-10]. The ground state 
of a clean 2D lattice with nested Fermi surface becomes 
unstable with respect to formation of antiferromagnetic 
spin gap under an arbitrarily small Coulomb interaction 
[11], and developes a Mott gap at strong interactions. 
Our recent studies of the DoS and of conductivity of 
a 2D electron gas on a lattice with substitutional im- 
purities [12,13] have shown that the commensurability 
effects at half-filling for the noninteracting case are op- 
posite to those obtained for ID systems, [14-16], so that 
the impurity scatterings with coherent Bragg reflections 
have been found to lower the electronic density of states 



around the Fermi level. Interaction effects in commen- 
surate weakly disordered 2D electronic systems have not 
yet been studied properly. Notice that, the recently ob- 
served metal-insulator transition (MIT) in 2D electronic 
systemSjWhich occurs at low temperatures (~ T < 2K), 
[17,18], cannot be explained in the framework of the 'con- 
ventional' localization theory and still is one of the puz- 
zling issues of central importance in the physics of dis- 
ordered systems. Measurements of resistivity in high- 
mobility Si-MOSFET's show that the insulator behav- 
ior at low particle densities, n < n c , crosses over to the 
metallic one as the particle density n reaches the crit- 
ical value n c — 9, 02 x 10 10 cm~ 2 ( for n > n c ). The 
observation of the MIT in GaAs/AlGaAs, [19], where e- 
e interactions are estimated to be weak, shows that the 
effect of correlations, although crucial, is not the only 
factor leading to MIT. 

In this Letter we study the effect of short-range repul- 
sive interactions on the one-particle DoS of a 2D square 
lattice with substitutional impurities for half-filled en- 
ergy band. We will show that a class of quantum correc- 
tions to the DoS, negative far from half-filling, changes 
its sign as the center of the band is approached. Such a 
behavior of the DoS is similar to that observed for the 
conductivity cr(T) [17-19]. A recent computation of the 
conductivity in the half-filled Hubbard model with dis- 
order also displays a change in the sign of da/dT as the 
system acquires particle-hole symmetry, [20]. 

The Hamiltonian of interacting electrons in the random 
field of substitutional impurities can be written in the 
Bloch-state representation in the following form, 
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e(p) = t{2 - cos(p x a) - cos(p y a)], 



(2) 



with t and a being the tunnelling integral for nearest- 
neighbor sites and the lattice spacing, respectively. 
Vi mp (q) and U(q) in Eq.(l) are the Fourier-transforms of 
a single impurity potential and of the short-range e-e in- 
teraction potential respectively; G is a reciprocal lattice- 
vector. pim P {q) = L~ 2 J2 a exp(iqR a ), where L is a lin- 
ear dimension of the system, and R Q is the coordinate 
of an impurity, randomly located on a lattice site. The 
impurity concentration is assumed to be small, and scat- 
terings on the ^-correlated impurity potential can be esti- 
mated in the framework of the Born approximation, [21]. 
For a metallic system ppl 3> 1 and crossed impurity lines 
of higher order in (ppl)' 1 are ignored; pp and I are the 
Fermi momentum and mean free path, respectively. 

By varying the band filling, Bragg reflections of 
the electronic wave off the Brillouin zone boundary 
are intensified for commensurate values of the electron 
wavelength, A, and the lattice constant, a. We will study 
the effects of Bragg reflections at half-filling, since they 
become essential as the middle of the band is approached. 
The Fermi surface for the energy dispersion given by 
Eq.(2) is flat at half filling, and its whole section is 
nested, with vectors Q = {±f , f }■ The perfect nesting 
of the Fermi surface gives rise to the following electron- 
hole symmetry relation for the electron dispersion with 
respect to Q: 



e(p + Q) - e F = -[e(p) - e F ], 



(3) 



where ep is the Fermi energy and ep = 2t for the half- 
filling case. Notice that introducing the next nearest- 
neighbor hopping term with amplitude t' destroys the 
perfect nesting of the Fermi surface. However, an optimal 
nesting takes place with Q* = 0.91Q for t'/t = 0.165, 
[22]. 

For a small band filling and far from the rational 
points, the Fermi surface of the model we are studying 
looks like a sphere, so that the effects of periodicity can be 
incorporated into the electronic effective mass. This situ- 
ation does not differ from that of an electron gas moving 
in the random field of impurities. The quantum correc- 
tions to the DoS for this model have been calculated in 
[8-10] . The lowest order quantum correction to the DoS 
resulting from the interactions in the diffusion channel 
can be expressed in the following form, [8,9]: 



Wr. t\ - J2d) [Ef(0,0)-2t7(p-p',0)] 



W(6,T) = -# 
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where D is the diffusion coefficient, D = VF 2 T ° , with 
vp and t being the Fermi velocity and the impurity 



relaxation time for Normal scattering, respectively, and 
h\{epmin{T , A }) . The bar over the in- 
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teraction potential U denotes the average over the Fermi 
surface, and the additional prefactor 2 in the Hartree 
correction term comes from the spin degeneracy. Inter- 
actions in the Cooper channel give the following contri- 
bution to the DoS far from half-filling, 
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where T c — e F exp(l/A) and A = p^ d ^U is the dimen- 
sionless interaction constant. Notice that the dynami- 
cal screening of a short-range interaction in the diffusion 
channel far from half-filling does not change the bare 
interaction potential considerably, whereas a short range 
potential is strongly renormalizcd in the Cooper channel, 
[21,24] as seen performing the summation of a ladder se- 
ries in the bare interaction, which replaces the logarith- 
mic energy or temperature dependence of the quantum 
correction by the rather weak double logarithmic depen- 
dence given by Eq.(5). 

A new kind of quantum corrections to the DoS takes 
place in the correlated disordered system due to the 
Bragg reflection, which enhances the number of electronic 
states at half-filling. 

As it is well known, the quantum interference correc- 
tions to the thermodynamic and kinetic coefficients come 
from the singular impurity ladder series referred to as 
the diffuson and the Cooperon blocks, [9]. The diffuson 
(Cooperon) block has the diffusion pole when the differ- 
ence q of the momenta of the electron and the hole (total 
sum k of the momenta of the electrons) and their ener- 
gies difference |w m | are small, i.e. the block acquires the 
pole in the diffusion regime when ql <C 1 (kl <C 1) and 

|^m|To < 1- 

New singular impurity blocks take place at half-filling 
with particle-hole symmetry, which are referred to as n- 
Diffuson and 7r-Cooperon, [12,13]. The 7r-Diffuson (n- 
Cooperon) has a diffusion pole at large oc Q momenta 
differences (total momenta) and small total energies of 
the electron and the hole (of two electrons). The dia- 
gram equations for the 7r-Cooperon, CV(q, iLu m ), and for 
the 7r-Diffuson, D 7T (q,iuj m ), are given in Fig.le,f. The 
bare Green's function Go is indicated in the diagrams by 
a straight line, and Gn(u, ie n ) = -. r-r-. — h-, — 5 — : 

uvr ' ie n — (e(p) — e f)+ 2^T signn n 

. The Green's function of an electron with large mo- 
mentum is represented by a dashed line. Notice that 
the dashing of the line has a physical meaning only for 
an impurity vertex or an interaction one when a large 
momentum transfer (oc Q due to the Bragg reflection) is 
involved in the scattering process, [12,13]. Therefore, the 
new selection of diagrams to be included in the summa- 
tion is performed according to the rule that a straight line 
in each impurity or interaction point joins a dashed line 
and vice versa. This rule is consistent with the fact that 
each scattering with large momentum transfer implies 
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also coherent Bragg reflection of the scattered electron 
on the boundary of the Brillouin zone. It is easy to see 
that the momentum conservation for impurity vertices 
with large momentum transfer is violated, so that they 
correspond to Umklapp scattering. As far as scatterings 
on point- like impurities are considered, these vertices 
are also characterised by t q . By summing the ladder se- 
ries in Fig. If, treating the 7r-scattering of an electron on 
impurities perturbatively, the following expression for the 
7r-Cooperon C ff (q, iu m ) is obtained: 



+ 



1 



-e„(w m - e„)) + 



2nr oP l d 

8(e n (uj m - e„)) 
(l + |a, m |T o ) 2 + (g0 2 -l 



(6) 



The expression for the 7r-Diffuson D v (q,iu) m ) is also 
given by Eq.(6) with the exception that q in the equation 
for D 7r (q, iu m ) denotes the differences of the momenta of 
a particle-hole pair. 
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FIG. 1. Bragg reflection induced contributions to the 
DoS at half-filling from (o), (o'); (c), (c') exchange and 
(b),(b');(d),(d') direct interactions. The diagrams obtained 
from the exchange ones by interchanging the straight and 
dashed lines joined in one interaction vertex with appropri- 
ate screening of the Coulomb interaction in particle - hole 
channel also give a contribution to the DOS. (e) and (f) are 
the impurity ladder series for the 7r-Diffuson,D^(q, iu> m ), and 
7r- Cooperon, Cr(q, iu> m ), respectively, (g) Diagram equation 
for the block D n (or CV), obtained by adding one impurity 
line to the 7r-Diflhson (or n- Cooperon). Heavy wave line and 
dotted line with cross denote Coulomb potential and impurity 
scattering, respectively. 

The additional contributions to the DoS in the middle 
of the band come from the diagrams shown in Fig.la- 
d'. The correction to the DoS, 5p(e,T), due to the e-e 
interactions is given by the following expression, 

2 f d 2 p 

6p(e,T) = --{Im J ^G'(p, ic„)S ee (p, "n)} ie „^ e 

(7) 

where S ee is the self-energy for first order quantum cor- 
rections to the DoS given by Figs.l(a)-(<f ). The expres- 
sion corresponding to the diagram Fig. la can be reduced 
to the following form after simple calculations, 

5p a (e,T) = 2rlpl d Im J £° ^f/ (q - 2 Q, 2e u ) 

(8) 



tanh ^ + tanh ^ 



[(i - ilot )2 + (qiy - iy 

The bare DoS p^ which is involved in the new cor- 
rection (8) at half-filling is again given by p^ d ^ = 
^yi t \n(e F min{T , l/|e|}), where the van Hove singular- 
ity is cut-off due to the i/2r term in the bare Green's 
function, [23]. Calculations of other diagrams in Fig.l 
are similar to the calculation of 5p a (e, T). Also, the con- 
tributions from the diagrams a'), b'), c') and d') in Fig.l 
are equal to those coming from a), 6), c) and d), respec- 
tively. Screening of the potential U in the diffusion and 
cooperon channels is realized according to the diagram- 
matic equations given in Fig. 2. 
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FIG. 2. Screening of the Coulomb interactions for elec- 
trons on a 2D half-filled square lattice in (a) particle-hole 
and (b)particle-particle channels, respectively. 
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In the middle of the band, it can be shown that the 
screening of a short-range interaction arises in the diffu- 
sion channel, instead of the screening of the potential in 
the Cooper channel which takes place in the electron gas 
model, [21,24]. It follows that a dynamical screening of 
the interaction in the diagrams (c) — (d') of Fig.l does 
not essentially change the value of the bare potential, 
and the contribution to the DoS from these diagrams is 
calculated according to Eq.(8), 



*JO U T , _ J2d) [£/(2Q,0) + 2t/(p-p' + 2Q,0)] 



m( 



1 



T max{e, T}' 



(9) 



The screening of the interaction in the diffusion channel is 
realized according to the diagrammatic equation in Fig. 2, 



pi 2d) [ln# -^(Ja^+i) + ^(i)] 

(10) 

The correction to the DoS from the diagrams in 
Fig.l(a, a') is calculated by putting Eq.(10) into Eq.(8). 
The contribution from these diagrams is given by the 
corresponding expression for 5p^\e,T), expressed by 
Eq.(5). The potential of the other diagrams in the dif- 
fusion channel, given by Fig. 1(6, b'), carries zero energy 
and large momentum ~ p' p" + 2Q. Therefore screen- 
ing is not effective and these diagrams give logarithmic 
contribution to the DoS. Therefore, the total quantum 
correction to the DoS in the diffusion channel can be 
presented as: 



6pi D He,T) = - 



1 



2ir 2 %D 



In 



In T c t /Ti 



p' + 2Q,0) 



2ir 2 HD 



ln( 



In T c /max{e,T} 
1 



+ 



T max{e,T} 



) (11) 



It is seen from Eqs.(9) and (11) that the Bragg reflection 
contribution to the DoS Sp v — Sp^ + Spi^ increases as 
e approaches the Fermi level. 

The total quantum correction to the DoS Sp(e, T) in 
the middle of the band can be expressed as Sp(e, T) = 
SpN(i,T)+SpT T (e, T). Both contributions to the DoS have 
the same energy or temperature dependences. However, 
they differ in sign and by the values of the interaction 
potentials U(0,0) and C/(2Q,0) . Since the short-range 
interaction is a screened Coulomb interaction, the po- 
tentials U(Q, 0) and f7(2Q,0) being included in Eq.(4) 
and Eq.(9) respectively cannot differ strongly from each 
other if the Thomas-Fermi screening number is of the 
order of the reciprocal lattice vector 2Q. As a result, the 
total quantum corrections to the DoS become positive 
at half-filling. Far from half-filling the Bragg reflections 
are destroyed and the dominant contribution to the DoS 



is the logarithmically decreasing Altshuler-Aronov con- 
tribution [8,9] given by Eq.(4). Notice that an anomaly 
in the one-particle DoS is experimentally observable by 
measuring the tunneling conductivity in a contact as a 
function of the bias voltage. This dependence completely 
reflects an energy dependence of the DoS. 

In conclusion, we showed that the DoS of a 2D dis- 
ordered lattice with flat Fermi surface is enhanced close 
to half-filling due to quantum corrections mediated by 
Bragg reflections. Since similar quantum corrections ex- 
ist for all commensurate points, the DoS seems to show 
oscillating behavior with maximum amplitude at half- 
filling. It is appropriate to emphasize that, although the 
Bragg reflection results in the decreasing the DoS of a 
non-interacting system, [12,13], it gives rise to an in- 
crease of the DoS in the presence of e — e correlations. 

The work was partially supported by ONR under 
Grant N00014-01-1-0427. 
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